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Abstract: The criteria used by International Union for Conservation of Nature (IUCN) for its Red List 
of Ecosystems (RLE) are the global standards for ecosystem-level risk assessment, and they have been 
increasingly used for biodiversity conservation. The changed distribution area of an ecosystem is one of 
the key criteria in such assessments. The S#pa bungeana grassland is one of the most widely distributed 
grasslands in the warm-temperate semi-arid regions of China. However, the total distribution area of this 
grassland was noted to have shrunk and become fragmented because of its conversion to cropland and 
grazing-induced degradation. Following the IUCN-RLE standards, here we analyzed changes in the 
geographical distribution of this degraded grassland, to evaluate its degradation and risk of collapse. Past 
(1950-1980) distribution areas were extracted from the Vegetation Map of China (1:1,000,000). Present 
realizable distribution areas were equated to these past areas minus any habitat area losses. We then 
predicted the grassland's present and future (under the Representative Concentration Pathway 8.5 
scenario) potential distribution areas using maximum entropy algorithm (MaxEnt), based on field survey 
data and nine environmental layers. Our results showed that the $. bungeana grassland was mainly 
distributed in the Loess Plateau, Hexi Corridor, and low altitudes of the Qilian Mountains and Longshou 
Mountain. This ecosystem occurred mainly on loess soils, kastanozems, steppe aeolian soils and 
sierozems. Thermal and edaphic factors were the most important factors limiting the distribution of S. 
bungeana grassland across China. Since 56.1% of its past distribution area (4.9104 km?) disappeared in the 
last 50 a, the present realizable distribution area only amounts to 2.2104 km*. But only 15.7% of its 
present potential distribution area (14.0*10* km?) is actually occupied by the $. bungeana grassland. The 
future potential distribution of S. bungeana grassland was predicted to shift towards northwest, and the 
total area of this ecosystem will shrink by 12.4% over the next 50 a under the most pessimistic climate 
change scenario. Accordingly, following the IUCN-RLE criteria, we deemed the $. bungeana grassland 
ecosystem in China to be endangered (EN). Revegetation projects and the establishment of protected 
areas ate recommended as effective ways to avert this looming crisis. This empirical modeling study 
provides an example of how IUCN-RLE categories and criteria may be valuably used for ecosystem 
assessments in China and abroad. 


Keywords: climate change; limiting factors; maximum entropy algorithm; potential distribution; realizable 
distribution; Red List of Ecosystems; Loess Plateau 


“Corresponding author: GUO Ke (E-mail: guoke @ibcas.ac.cn) 

Received 2018-06-13; revised 2019-09-13; accepted 2019-10-10 

© Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Science Press and Springer-Verlag GmbH Germany, 
part of Springer Nature 2020 


http://jal.xjegi.com; www.springer.com/40333 


JOURNAL OF ARID LAND 


1 Introduction 


The massive loss of ecosystem services induced by inappropriate human activities in the last few 
decades has undermined the welfare of human beings (Costanza et al., 2014), with 37% of the 
global ecoregions since converted to human-dominated uses (Hoekstra et al., 2005). Considering 
their variations in spatial range, species composition, and community structure and functioning 
under environmental change and human interference, many ecosystems have become endangered 
or have collapsed already (Miles et al., 2006; Burns et al., 2015). Adopting a conservation 
outlook and strategy, priorities should be first set by identifying those ecosystems currently most 
at risk. To this end, the International Union for the Conservation of Nature (IUCN) has 
developed a set of categories and criteria for assessing the risk posed to ecosystems and for the 
efficient conservation of these ecosystems, known as the Red List of Ecosystems (RLE) (Bland 
et al., 2015). This assessment at the ecosystem level provides a useful body of information to 
support the biodiversity conservation, macroeconomic planning and management decisions 
(Keith et al., 2015). Importantly, the categories and criteria of an IUCN-RLE assessment can be 
applied to any habitat type, biotope, ecosystem, vegetation type and ecological community; all of 
them are considered suitable assessment units and regarded as operational synonyms of the 
ecosystem type (Bland et al., 2015). Since their introduction, the IUCN-RLE categories and 
criteria have been widely applied across a variety of ecosystem types and geographical areas 
(Bland et al., 2015; Burns et al., 2015; Ma, 2017). 

China's natural ecosystems have experienced long-term anthropogenic activities, with many of 
its ecosystem types collapsed and many others only persisted in a highly fragmented state. The 
Ministry of Ecology and Environment of China has set forth a series of projects to compile a 
national red list of ecosystems for conserving their immense biodiversity and sustaining 
ecosystem services. The vast grassland of northern China is one area of focus in this assessment. 
In particular, Stipa bungeana grassland is the most representative plant community of the 
warm-temperate grasslands that span the Palearctic biome, with its major distribution area found 
on the Loess Plateau of China (Chen et al., 2014; Cheng et al., 2015; Wesche et al., 2016; Guo et 
al., 2018). Due to large-scale grassland reclamation (conversion to arable land) and degradation 
caused by over grazing, S. bungeana grassland has rapidly declined in distribution and become 
increasingly fragmented over the past century (ECVC, 1980; ECVIM, 1985; Jiao et al., 2008). 
Meanwhile, human-induced climate change, including rising temperatures and changed 
precipitation patterns (IPCC, 2013), is expected to drive substantial shifts in species’ ranges and 
to alter the distribution pattern of communities (Parmesan and Yohe, 2003; Root et al., 2003), 
perhaps even causing extinctions of species and ecosystems (Walther et al., 2002). S. bungeana 
plants are sensitive to high temperatures and elevated atmospheric CO2 concentrations (Song et 
al., 2016; Wang et al., 2017; Shi et al., 2018). Hence, global warming may further reduce the 
geographical distribution of S. bungeana grassland and increase its risk of collapse, by decreasing 
niche diversity and increasing its susceptibility to competition and stochastic perturbations 
(Harpole and Tilman, 2007; Keith et al., 2013; Bland et al., 2015). Thus, rigorously assessing the 
collapse risk of S. bungeana grassland is urgently needed to inform and support the 
decision-making process for the conservation and sustainable management of grasslands in 
China. Doing this also would exemplify how IUCN-RLE categories and criteria could be valuably 
used for conducting rigorous ecosystem assessments in China. 

Changes in distribution area constitute a key criterion in the IUCN-RLE assessment 
framework. Accordingly, in this study, we assessed the status of S. bungeana grassland in China 
by focusing on its changed distribution areas. We first addressed the following two questions: (1) 
how much of the original S. bungeana grassland habitat was lost during the last few decades and 
how would its distribution range shift under future climatic changes? and (2) what is the gap 
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between the potential distribution area of the grassland and its realizable one and which factors 
are most limiting this distribution? Then, considering its distribution area changes found and 
other IUCN-RLE criteria, we assigned the S. bungeana grassland to an IUCN-RLE, and 
recommended conservation measures for this ecosystem in semi-arid regions of China. 


2 Materials and methods 


2.1 Description of S. bungeana grassland and criteria used to assess its collapse risk 


S. bungeana is a xerophytic and thermophilic bunchgrass plant distributed in eastern Eurasian 
grasslands and forest steppe areas. When this plant species is locally dominant, it forms the S. 
bungeana grassland. In China, this grassland occurs mainly on the Loess Plateau, at altitudes 
ranging from 631 to 3174 m, where the mean annual precipitation is 374 mm and the annual mean 
temperature is 6.7°C. This grassland has a relatively simple species composition (richness of 9 
species/m7): Cleistogenes squarrosa, Agropyron cristatum, Heteropappus altaicus and Astragalus 
scaberrimus being the most common species in the community. Much of this grassland has been 
converted into arable lands for crops, limiting its current distribution to steep slopes and 
abandoned croplands. 

Two criteria (A and B) were used to assess the collapse risk of China's S. bungeana grassland: 
the reduction of its natural distribution and the extant restriction to its distribution. To determine 
the former, the distribution areas of S. bungeana grassland were compared over a time series 
(past, present and future), while the latter was based on its extent of occurrence (EOO) in the 
presently distributed area, as well as its area of occupancy (AOO) (Fig. 1). 
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Fig. 1 Framework used to assess the collapse risk of Stipa bungeana grassland based on changes in its spatial 
distribution. Vegetation Map of China (VMC; 1:1,000,000) was used to extract the pertinent vegetation 


information for the period 1950-1980. This map was originally generated from field survey data for this period. 
ENMs is the abbreviation for ecological niche models. 
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2.2 Past (1950-1980) and present (2010) realizable distributional data 


The "past" refers to a timeframe encompassing the last 50 a to date in the IUCN-RLE assessment 
(Keith et al., 2015). Here, the past distribution data were derived from the Vegetation Map of 
China (VMC; 1:1,000,000) (ECVMC, 2007), which is drawn from extensive vegetation survey 
data collected in 1950-1980. The land area corresponding to the S. bungeana grassland was 
extracted in ArcGIS v10.0 (ESRI Inc., Redland, CA, USA). 

The habitat losses of S. bungeana grassland were removed from the past realizable distribution 
to acquire its present distribution. Land use and land cover (LUCC) data for 2010, provided by 
Data Center for Resources & Environmental Sciences, Chinese Academy of Sciences, were used 
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to examine the land type covered by S. bungeana grassland from the VMC. Six land cover types 
were identified, among which "grassland" was considered the only suitable habitat for the S. 
bungeana grassland, whereas the other five were deemed unsuitable (i.e., cultivated land, 
woodland, water bodies, building land and unutilized land (sand or bare land)). Importantly, the 
distribution of S. bungeana grassland can be reduced by not only losses to the five unsuitable land 
covers, but also to grassland that becomes dominated by other species such as Artemisia wellbyi. 
The former outcome was generated with LUCC data and the latter determined via field surveys in 
the recorded distribution areas. 


2.3 Predicting the potential distribution 


2.3.1 Model selection 


Ecological niche models (ENMs), also known as species distribution models (SDMs), are being 
increasingly used to guide effective conservation policy and action (Guisan et al., 2013). Now 
widely used to estimate threshold ecological requirements, they can map the potentially suitable 
habitat of both species and communities (Franklin, 2009; Peterson et al., 2011; Phillips and Elith, 
2013; Prieto-Torres et al., 2016). ENMs rely on Hutchinson's super-volume theory (Hutchinson, 
1957) and emphasize the role of abiotic factors for governing distribution of flora and fauna 
(Guisan and Thuiller, 2005). 

Dozens of such models exist, with MaxEnt (maximum entropy algorithm) currently offering 
perhaps the best modeling approach (Moisen and Frescino, 2002; Segurado and Araujo, 2004; 
Elith et al., 2006). This model can make use of presence-only spatial data, as well as both 
categorical and continuous environmental variables (Phillips et al., 2004, 2006). Specifically, 
MaxEnt uses a maximum -likelihood method to generate the probability distribution in space over 
the pixels. This probability is termed "gain", which is defined as the average log probability of the 
presence samples minus a constant that makes the uniform distribution have zero gain. The gain 
starts at zero and increases as the algorithm enhances the probability with more sampled 
locations. Through an iterative procedure, the gain would continue to increase until it fell below 
the convergence threshold. The gain thus represents the likelihood of samples occurring. For 
example, if the gain in a pixel is 3, the likelihood of distribution in this pixel is exp (3), which is 
20.1-fold that of a random background pixel (Yost et al., 2008). In this way, MaxEnt produces a 
continuous prediction, with values ranging from 0 to 1: a cell scored 1 is the most suitable habitat 
for distribution, while a cell scored 0 is the least suitable. 

Elith et al. (2010) had used exploratory techniques to map the abiotic factors limiting the 
distribution of a given species. The smallest value in each grid cell was chosen as the limiting 
factor after calculating a multivariate environmental similarity surface (MESS) (Elith et al., 
2010). That analysis underscored the most important ecological factor drive the obtained 
predictions. The MESS value of a point P is calculated as follows: 


(p, —min,)/(max,—min,)x100 f, =0 

x f, 0< f, <50% 
2x (100- f,) 50% < f, <100% |’ 
(max, — p,)/(max, —min,)x100 f, =100% 


MESS = min(MESS, ) = min (1) 


where i is the i predictor variable; min; and max; are the minimum and maximum values of 
variable i over the reference point set, respectively; p; is the value of variable i at point P; and fi 
(%) is the percentage of reference points whose value for variable i is smaller than pj. 

2.3.2 Plot presence data 

We conducted extensive field work during the summer (July to September) from 2011 to 2017 on 
China's steppe biomes of the Northeast China Plain, the Inner Mongolian Plateau, the Loess 
Plateau, the Tibetan Plateau, and the mountainous areas of Northwest China. A total of 899 plots 
(20 mx50 m) with 3-5 quadrats (1 mx1 m) in each plot were set up. For each plot, its GPS 
coordinates, geomorphic type, slope aspect and inclination, land use information and soil type 
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were recorded. All living vascular plants and dead annuals within the quadrats were recorded at 
the species level by their rooted presence. Then their individual heights (cm) were measured on a 
per species basis, from which coverage by each plant species was estimated as a percentage of the 
quadrat area. For each species in a quadrat, we also clipped its aboveground parts at the soil 
surface, then oven-dried these samples at 65°C to a constant weight to obtain their respective final 
dry biomass (g/m*). The importance value (IV) of each species (Curtis and McIntosh, 1951) was 
calculated this way: IV=(relative height+relative cover+relative biomass)/3, and was used to 
quantify their rank dominance in the vegetation. The highest IV was deemed the dominant species 
in the plot and used to name this community. In this way, 103 sites were found dominated by S. 
bungeana in our survey work. Along with an additional eight sites reported by Cheng et al. 
(2010), the locations of these 111 vegetation sites of S. bungeana grassland were used to predict 
its distribution with the MaxEnt (Fig. 2). 
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Fig. 2 Locations of the 111 vegetation sites dominated by S. bungeana grassland that were used for predicting 
its distribution in MaxEnt (a), the current status of S. bungeana grassland in abandoned cropland (b) and artificial 
Caragana microphylla shrubland where S. bungeana is growing between shrubs (c) 


2.3.3 Environmental variables 


Both climatic and edaphic variables were used to discern the ecological niche of the studied 
community. First, a set of 19 climatic variables (at 30 s resolution, approximately 0.9 kmx0.9 km) 
from WorldClim Version 2 was downloaded, representing the averages for the 1970-2000 period 
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(Fick and Hijmans, 2017). Strong collinearity of among these variables in MaxEnt can influence 
the final modeling result (Phillips et al., 2017), so a Pearson's correlation analysis was applied to 
them to reduce any redundant information in the data. Specifically, only one variable was retained 
if the correlation coefficient of pair-wise variables was >0.8. This led to 8 of the 19 climatic 
variables chosen to represent the thermal and moisture factors, i.e., mean diurnal range (MDR), 
isothermality (ISO), temperature seasonality (TS), maximum temperature of the warmest month 
(MTWM), minimum temperature of the coldest month (MTCM), precipitation of the wettest 
month (PWM), precipitation of the driest month (PDM) and precipitation seasonality (PS). The 
sole edaphic variable used here was soil type, taken from Soil map of China (at 1.0-km 
resolution) based on a generic classification of soil. The Soil map of China was downloaded from 
the Soil Data Center, National Earth System Science Data Sharing Infrastructure, National 
Science & Technology Infrastructure of China (http://soil.geodata.cn). To ensure that all climatic 
and edaphic variables shared the same geographical bounds and cell size, we overlaid them in 
ArcGIS v10.0 to create a mutual area with a cell size of 42 s (1.3 kmx1.3 km). Hence, in this way, 
China's land map was divided into a grid comprising 5806 columns and 3290 rows. 


2.3.4 Parameters’ run in MaxEnt and model evaluation 


With the S. bungeana presence data and nine environmental layers (eight climatic layers plus one 
soil type layer) prepared, models were generated in MaxEnt v3.4.0 by randomly sampling 75% of 
the locality records as training data, leaving the remaining records to serve as testing data (for 
model validation). We also ran 1000 iterations with no extrapolations to avoid an artificial 
extrapolation from extreme values of the ecological variables (Elith et al., 2011). The cloglog 
output format was used, as recommended by Phillips (2017), to express the continuous presence 
probability on a scale from 0 to 1. All other parameters were set to their default settings. 

The area under the curve (AUC) of the receiver operating characteristics curve was relied upon 
to quantify the performance of the MaxEnt (Anderson et al., 2003; Elith et al., 2006; Elith et al., 
2011). An AUC value ranges from 0.0 to 1.0, according to which a model's performance can be 
roughly classified as "poor" (0.5—0.6), "fair" (0.6—0.7), "good" (0.7—0.8), "very good" (0.8—0.9) 
and "excellent" (0.9-1.0). The relative importance of each variable was determined by its 
contribution to the model and a jackknife test of the regularized training gain (Phillips, 2017). 


2.3.5 Suitable categories 


Normally, the threshold used to shift a model's form, from continuous probability to binary, is set 
by the "maximum sensitivity plus specificity" value provided by MaxEnt (Brites-Neto and 
Duarte, 2015). This particular value for the S. bungeana grassland ecosystem was 0.3. The 
ensuing probability values on the distribution map were expanded 100 times to facilitate several 
interpretations. We then divided the map into four suitable categories based on these latter 
probability values: core area (75-100), medium area (50-75), marginal area (30-50) and 
unsuitable area (0—30). To bolster confidence in our predictions, we considered the core and 
medium areas as distribution areas of S. bungeana grassland, but the marginal and unsuitable 
areas as its absence area. To eliminate further unsuitable habitat areas, the current potential 
distribution map was filtered by LUCC maps. 


2.3.6 Future climate change scenario and predicting the future distribution 


With the goal of predicting future trends for S. bungeana grassland's distribution in China, we 
predicted its distribution in 2070 (ie., averaged over years for 2061-2080) based on the 
Representative Concentration Pathway (RCP) 8.5 scenario. The RCP 8.5 reflects relatively slow 
income growth and modest improvements in technology and energy-use intensity, leading to 
higher demand for fossil fuels and rising greenhouse gas emissions in the long term in the 
absence of enacted climate change policies (Riahi et al., 2011; IPCC, 2013). Compared with other 
RCPs, the RCP 8.5 is the worst-case scenario (highest greenhouse gas emissions), so the results 
under this scenario could be interpreted as pessimistic. Under the RCP 8.5, the global atmospheric 
CO: concentration is projected to reach 677 cm?/m? by 2070 (IPCC, 2013). 

Future climate layers were generated using global climate models (GCMs), many of which are 
currently available. To make inferences that are robust under a wide range of possible futures, it is 


QIAO Xianguo et al.: Assessing the collapse risk of Stipa bungeana grassland in China based on... 


advised to use an ensemble of these models, to smooth out their variability and to provide a 
consensus result, instead of choosing one particular model or pathway (Moss et al., 2010; IPCC, 
2013). Here, five GCMs from five modeling centers in four countries formed this ensemble: 
BCC-CSM1-1 (Beijing Climate Center Climate System Model v1, China), CCSM4 (Community 
Climate System Model v4, the National Center for Atmospheric Research, USA), GISS-E2-R 
(Goddard Institute for Space Studies Model E2 Earth System Model with the Russell ocean 
model, USA), HadGEM2-AO (Hadley Centre Global Environment Model v2, Met Office Hadley 
Centre, UK), and MIROC-ESM-CHEM (Atmospheric Chemistry Coupled Version of Model for 
Interdisciplinary Research on Climate-Earth System, Japan). The bioclimatic data used by all of 
them were downloaded from Worldclim (http://www.worldclim.org/). 

Five future distribution models were generated in MaxEnt by the five GCMs. The parameters 
and standards for category division were the same as those used for predicting the present 
potential distribution of S. bungeana grassland. The final, future distribution model was 
integrated by averaging the five models generated by MaxEnt. Lacking a land cover map of the 
future, the resulting future distribution cannot be filtered by land cover or land use data; hence, 
the distribution change from present to future was analyzed by using maps without filtration. The 
"unchanged area" thus refers to the same area found suitable in the present and future; the "gained 
area" was calculated after removing the unchanged area from the future distribution area; and the 
"lost area" was calculated by deducting the unchanged area from the present distribution area. 


2.4 Formally assessing the risk of ecosystem collapse 


Eight categories and five criteria are currently set by the IUCN-RLE, for assessing the risk of an 
ecosystem collapse. These are collapsed (CO), critically endangered (CR), endangered (EN), 
vulnerable (VU), near threatened (NT), least concern (LC), data deficient (DD) and not evaluated 
(NE). They are defined by five quantitative and qualitative criteria: declining distribution (A), 
restricted distribution (B), environmental degradation (C), disruption of biotic processes and 
interactions (D), and quantitative risk analysis (E). The overall ranking of an ecosystem is 
determined by the most severe ranking of these five criteria, of which two (A and B) are directly 
related to spatial distribution, each having three sub-criteria. So, in this study, the S. bungeana 
grassland was assessed by criteria A and B. 

2.4.1 Criterion A 


This criterion requires an assessment of decline in geographical distribution across three time 
periods: past (over the past 50 a), future (over the next 50 a) and historical (since 1750) 
reductions. Past and future declines share the same thresholds: a decline >80% denotes CR, 250% 
corresponds to EN and >30% indicates VU; but if neither of these is met, then a category of LC 
(Least Concern) applies. Historical decline could not be assessed in this study, however, because 
the required distribution data for S. bungeana grassland remain unavailable. 


2.4.2 Criterion B 


When an ecosystem is strongly restricted, its risk of collapse is higher when confronted with a 
spatially explicit threat or catastrophe (Rodríguez et al., 2015). So, this criterion serves to identify 
those ecosystems whose distribution is so restricted they could collapse from a stochastic 
occurrence of a few interacting threatening events. To fully apply criterion B, assessors must 
calculate the EOO, AOO and the number of locations from a suitable distribution map. EOO is 
the area contained within the smallest polygon encompassing all known, inferred or projected 
sites of an ecosystem's occurrence at present. An EOO area of <0.2x10* km? is CR, <2.0x10* km? 
is EN and <5.0x10* km? is VU. AOO is the area within EOO occupied by an ecosystem type, 
measured in 10 kmx10 km grid squares: a number of <2 grid squares is CR, <20 is EN and <50 is 
VU. The number of locations refers to a small number of these, generally <5, prone to the effects 
of human activities and stochastic events in an uncertain future, thus capable of collapse or 
becoming CR within a very short timeframe. This subcriterion can only lead to a VU category. 
Only the EOO and AOO were applied here because of insufficient data for such locations. 
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3 Results 


3.1 Past and present realizable distributions of S. bungeana grassland 


The past distribution area of S. bungeana grassland extracted from the VMC and representative 
for 1950-1980, totaled 4.9x10+ km?. By the year 2010, 61.2% (3.0x10* km?) of this area had 
persisted as natural grassland, with the rest mostly converted to cultivated land (27.2%) and 
woodland (7.0%). According to the VMC, one could find extensive S. bungeana grassland areas 
distributed in the basin along the Yarlung Zangbo River and its tributaries in the Tibetan Plateau 
before 1980. Surprisingly, when thorough surveys were conducted in those areas during summer 
time of 2013, 2018 and 2019, not a single vegetation patch dominated by S. bungeana could be 
found, just a few of its individuals; instead, now prevailing in those areas were A. wellbyi and 
Tripogon bromoides communities. Consequently, the present realizable distribution area of S. 
bungeana grassland was in actuality 2.2x10* km? (44.9% of the past distribution) after 
eliminating this large affected area in the Tibetan Plateau (Fig. 3a). 
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Fig. 3 The past and present realizable distributions (a), limiting factors of the distribution (b), present potential 
distribution (c) and future potential distribution under Representative Concentration Pathway (RCP) 8.5 (d) of 
Stipa bungeana grassland in China. ST, soil type; PWM, precipitation of the wettest month; PDM, precipitation of 
the driest month; PS, precipitation seasonality; MDR, mean diurnal range; ISO, isothermality; TS, temperature 
seasonality, MTWM, mean temperature of the warmest month; MTCM, mean temperature of the coldest month. 


3.2 Present potential distribution of S. bungeana grassland and its driving factors 
3.2.1 Model performance and limiting factors for the present potential distribution 
The applied model generated high values for the AUC test (training data: 0.984; test data: 0.977), 
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indicating its predictions were highly robust. Soil type (ST), PWM and MTCM were the most 
important factors, together accounting for 74.0% of the variation in the predicted distributions 
(Table 1). Factors with respective contributions of >10.0% could be categorized into thermal 
(MDR, MTCM and MTWM), moisture (PWM) and edaphic (ST) groups, with the thermal 
(39.4%) and edaphic (37.8%) factors clearly most prominent, whereas the moisture factor 
explained less of the distributional variation (22.8%). 

Table 1 MaxEnt model-predicted climatic ranges for the present suitable distribution area of Stipa bungeana 
grassland across China, and the respective contribution of nine environmental variables to these predictions 


Present potential suitable habitat category 


Variable Contribution (%) 


Core area Medium area Marginal area Unsuitable area 
PWM (mm) 27 45.0-135.0 41.0-140.0 32.0-159.0 3.0-1244.0 
PDM (mm) 0.0 0.0-6.0 0.0-18.0 0.0-18.0 0.0-206.0 
PS 0.1 78.3-121.1 37.7-122.5 35.2-141.5 19.1-155.7 
MDR (°C) 11.6 10.2-16.1 9.2-16.2 8.7-16.2 1.1-20.1 
ISO (%) 2.9 24.2-39.6 23.4-41.8 21.7-45.4 4.4-76.3 
TS 22 8.0-13.3 6.9-13.5 6.4-14.4 0.0-18.3 
MTWM (°C) 9.9 18.0-31.2 17.0-31.8 14.5-32.1 —15.7-40.3 
MTCM (°C) 12.8 -22.3- -9.2 -23.5- -7.2 —24.1--4.3 —48.1- -14.5 


ST 37.8 - - - - 


Note: PWM, precipitation of wettest month; PDM, precipitation of driest month; PS, precipitation seasonality, MDR, mean diurnal 
range; ISO, isothermality; TS, temperature seasonality, MTWM, maximum temperature of the warmest month; MTCM, minimum 
temperature of coldest month; ST, soil type. -, no value. 


PWM was by far the strongest factor limiting the distribution of S. bungeana grassland in 
China. Specifically, it explains why 65.9% of China's land area (in the east, south and northwest) 
lacks S. bungeana grassland. MTCM was the second most critical factor, limiting the grassland's 
distribution to the Tibetan Plateau and Northeast China, accounting for 26.2% of China's land 
area (Fig. 3b). 

3.2.2 Present potential distribution area and its climatic range 


The model generated a suitability map of S. bungeana grassland (Fig. 3c). According to the 
probability values obtained, four suitable categories (core, medium, marginal and unsuitable 
areas) were defined. The total distribution area of S. bungeana grassland (i.e., core area+medium 
area) was 28.0x10* km?, thus accounting for 2.9% of China's land area. This ecosystem is 
sustained mainly on the Loess Plateau, Hexi Corridor, and at low altitudes of the Qilian 
Mountains and Longshou Mountain, with very few patches of it scattered in the Gonghe Basin, 
south of the Chang Tang Plateau and Liaohe Plain. Overall, 86.9% of the found suitable areas 
occurred on the Loess Plateau, of which 13.5% and 22.9% there respectively served as core and 
medium areas. Not all of the realizable distribution was included in the present potential 
distribution, for which 59.2% and 81.8% of the past and present realizable distribution were 
matched to, respectively, with the rest of the past realizable distribution mainly distributed in the 
basin along the Yarlung Zangbo River, and the rest of the present realizable distribution mainly in 
the Northeast Plain. Only two field survey plots are beyond the east boundaries of the present 
potential distribution area. 

Ranges for the environmental factors of the distribution's core area are summarized in Table 1, 
for which the most common soil types were loess soils (46.0%), kastanozems (18.7%), sierozems 
(9.0%) and steppe aeolian soils (1.1%). Filtering this by land cover data left only 50.9% 
(14.0x10* km?) of the present potential distribution area available to sustain the S. bungeana 
grassland, or 1.5% of China's land area. 


3.3 Future potential distribution of S. bungeana grassland under climate change scenario 
(RCP 8.5) 


Changes to the future distribution area of China's S. bungeana grassland ranged from —16.4% to 
0.6%, with four of the five ensemble models predicting a decrease, whereas the other predicted an 
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increase; for 2070, their arithmetic mean was —12.4% (Fig. 3d). This predictive modeling exercise 
revealed that 70.9% of the present distribution area would remain unchanged, with gained and 
lost areas of 16.7% and 29.1%, respectively. Habitat losses occurred mainly on the eastern and 
southern borders of the grassland's present distribution area, which would cause it to narrow in 
range. The gained habitats would mostly occur in the Qaidam Basin and Tianshan Mountains, in 
addition to the west part of the Liaohe Plain; that is, in the western and northeastern borders of 
the present distribution areas. Evidently, the S. bungeana grassland distribution should shift 
northwest in the future, with a modest reduction in its total area. 


3.4 Collapse risk of S. bungeana grassland 


In just the past 50 a, S. bungeana grassland's total distribution area has declined by 65.1%, with a 
further 12.4% reduction driven by climate change expected over the next 50 a. Based on these 
changes in distribution area, and assuming no further changes to it induced by human activities 
for the predicted period (i.e., no further conversion of natural grasslands to croplands, or vice 
versa), we ranked the S. bungeana grassland as endangered (EN) according to criterion A of the 
IUCN-RLE. Given the patchy occurrence of S. bungeana grassland, 2.2x10* km? in 301 of the 10 
kmx10 km grid cells, this ecosystem can be ranked as vulnerable (VU) as per criterion B. Overall, 
the S. bungeana grassland in China is deemed as endangered (EN). 


4 Discussion 


4.1 Climatic factors shaping the distribution range of S. bungeana grassland 


Climatic factors are fundamental for shaping the distribution of vegetation (Pan et al., 2003). 
Either thermal (Kira, 1976; Fang and Lechowicz, 2006; Li et al., 2016a) or moisture (Stephenson, 
1990, 1998; Byrne et al., 2017) factors can play primary roles in controlling the boundaries of 
plant species or communities. S. bungeana plants inhabit warm-temperate zones, reaching the 
southern limit of steppes in China, and the most outstanding trait of this species is its heliophilic 
idiosyncrasy (ECVIM, 1985; Wu and Raven, 2013). The ranking of variables by MaxEnt 
confirmed that thermal factors play a more prominent role than moisture factors in the 
distribution of S. bungeana grassland, a result consistent with it having larger variation (i.e., a 
greater coefficient of variation (CV=standard deviation/mean)) in precipitation than temperature. 
Similar conclusions were reached for different plant species studied in China (e.g., Fang and 
Lechowicz, 2006; Li et al., 2016b). However, Wen et al. (2008) found that S. bungeana's 
distribution along the Yanhe River catchment was mainly controlled by annual evaporation and 
temperature seasonality. This discrepancy implies that predictions of S. bungeana's spatial 
distribution areas may vary with different climatic layers used and presence data covering 
different parts of its native range. Nevertheless, the S. bungeana grassland was most commonly 
found on the Loess Plateau (ECVC, 1980; Cheng et al., 2010), which suggests that soil type is a 
key contributing factor to consider when building its distribution model. We conclude that 
thermal and edaphic factors are crucial for shaping the distribution of S. bungeana grassland 
across China. 


4.2 Gap between realizable and potential distributions 


A variety of factors can make the realizable distribution range much smaller than the present 
potential one, these include habitat availability, dispersal limitation, interspecific competition, 
extinction events and historical constraints (Pulliam, 2000; Guisan and Thuiller, 2005). The ratio 
of the realizable range size to potential range size, i.e., range filling, is a measure useful for 
quantifying how well species fill their potential ranges (Gaston, 2003). For example, the low 
range filling (38.3%) of many tree species in Europe is considered driven primarily by 
geographical dispersal constraints on their post-glacial expansion and the climate (Svenning and 
Skov, 2004). In our study, upon excluding unsuitable land cover, range filling of S. bungeana 
grassland was estimated to be 15.7% (i.e., 2.2x10* km?/14.0x10* km?). Grassland degradation 
and relief variation could well be the most decisive factors for this extremely low range filling. 
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Overgrazing is a major human-imposed factor that initiates the retrogressive succession of 
grassland communities, which eventually leads to the turnover of species composition and 
grassland degradation (Mclendon and Redente, 1990; Young et al., 2001; Akiyama and 
Kawamura, 2007; Mirzabaev et al., 2016). Over time, as overgrazing continues unabated, the S. 
bungeana grassland may become a community instead dominated by Artemisia frigida or 
Potentilla acaulis plants. In recent decades, grassland degradation has been a common 
phenomenon in China, along with human population growth and its associated activities, such 
that a large proportion of plant species’ potential distribution area is now occupied by degraded 
community types (Fu, 1989; Akiyama and Kawamura, 2007). Differences in relief often bring 
with them differing dominant species. Several plant communities often coexist as mosaics within 
one landscape due to the complex landforms in the distribution area of S. bungeana grassland 
(Fu, 1989; Mclendon and Redente, 1990). For example, when S. bungeana grassland occurs on 
north-facing loess slopes, its south-facing slopes are usually covered by Stipa breviflora 
grassland, with the gully occupied by Artemisia giraldii grassland while Thymus mongolicus 
covered the stony surface of hills (ECVC, 1980). Evidently, small patches of realizable S. 
bungeana grassland can strongly or disproportionately influence the potential distribution map 
by making it more extensive. Conversely, S. bungeana grassland can also occur in places where 
the regional climate does not seem optimal, but the local climate, such as the sunny slope of a 
shrub prevailing landscape, might suit the S. bungeana grassland quite well. We surmise this 
could explain the presence of realizable distribution and field survey plots found outside the 
potential suitable area. 


4.3 Role of revegetation approaches in S. bungeana grassland's distribution 


Planting trees or grasses and abandoning cropland are common ways to achieve vegetation 
restoration, especially on the Loess Plateau. Yet the former approach (planting trees or grasses) 
does not seem to favor the recovery of S. bungeana grassland, but undoubtedly it has contributed 
immensely to increased vegetation productivity there (Feng et al., 2016). Shrub-like trees, such as 
Caragana microphylla and Hippophae rhamnoides, now occur in many parts of this region; 
however, their widespread planting has generated an unsustainable water supply and conditions 
that preclude the natural succession of grasslands (Chen et al., 2002). Alternatively, when 
cropland is abandoned for natural restoration, it must normally go through five or six stages until 
a vegetation climax state is reached. This process of secondary succession goes through 
communities dominated by Artemisia scoparia, Lespedeza davurica, Agropyron cristaum, L. 
davuricat+Artemisia sacrorum, Bothrichloa ischaemum, and finally S. bungeana (Hao et al., 
2005). But it usually takes 20 to 40 a before S. bungeana attains dominance in the community. 
Most of the abandoned croplands are in late successional stages, dominated by forbs or 
rhizomatous grasses. Therefore, it is plausible that much of that abandoned land will become S. 
bungeana grassland in the near future. 


4.4 Reasons for losses in the past distribution in the basin along the Yarlung Zangbo River 


It is quite perplexing that 1.0x10* km? of S. bungeana grassland in the Tibetan Plateau 
disappeared entirely in just the past 50 a (Zhao et al., 2018). Besides those areas occupied by 
reclamation and building land, 80% of this grassland was displaced by a shrub-grassland 
continuum dominated by A. wellbyi and T. bromoides. We propose this vast disappearance is the 
combined outcome of global climate change effects and natural vegetation succession. The 
Tibetan Plateau is one of the world's most sensitive areas for responding to global climate change 
(Liu et al., 2009). Already, it has experienced significant warming since the mid-1950s, and this 
warming tendency may continue under future climate change scenarios (Liu and Chen, 2000). 
Firstly, since the main factors limiting the distribution of S. bungeana's grassland in the basin 
along the Yarlung Zangbo River are thermal ones, an increasing temperature could have made the 
past suitable distribution area unsuitable. Secondly, the cyclic vegetation succession of sandy 
grassland includes two stages: a positive one and a negative one. The positive part starts with 
pioneer species colonizing drifting sand dunes up to Artemisia ordosicat+Artemisia 
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sphaerocephalai on semi-fixed sand dunes, reaching to A. ordosicat+S. bungeana+mosses on 
fixed sand dunes, and eventually the S. bungeana grassland. The negative part begins with 
desertification and reverts the positive succession backwards to A. sphaerocephalai-dominated 
communities or drifting sand dunes. This positive (or negative) phenomenon has been well 
studied in the Mu Us Sandy Land of the Mongolian Plateau (Guo, 2000). When vast portions of 
the Tibetan Plateau grassland are degraded because of overgrazing, desertification inevitably 
follows (Harris, 2010). Those plant species tolerant of sandy environments, such as A. wellbyi and 
T. bromoides, are, in turn, favored to invade and become dominant. In other words, it represents a 
retrogressive succession similar to the one in Mu Us Sandy Land, only with different endemic 
species involved. The elaborate features and driving forces involved in these grassland 
successional dynamics merit further study. 


4.5 Collapse risk and recommendations for the conservation of S. bungeana grassland 


The widely distributed S. bungeana grassland has been largely destroyed, converted either into 
cropland (Fu, 1989; Feng et al., 2016) or degraded into other grassland types (ECVC, 1980). 
Those vegetation communities still dominated by S. bungeana are highly fragmented, persisting 
as small patches on villages’ cemetery, on steep slopes or on abandoned lands in most of its 
distribution area on the Loess Plateau. The only one well preserved site with a continuous and 
extensive distribution of S. bungeana grassland ecosystem type lies protected in the Yunwu 
Mountain Natural Reserve in Ningxia (Nature Reserve Management Office of Yunwu Mountain 
in Ningxia, 2001). Thus, as an endangered plant community, the S. bungeana grassland calls for 
immediate on-the-ground action for its conservation. Sowing S. bungeana seeds in those areas on 
the Loess Plateau, particularly where the habitat suitability for this species is high both at present 
and in the future, is a strongly recommended measure for vegetation restoration, soil conservation 
and forage production, as S. bungeana tolerates grazing (Wang, 1989) and is a drought-resistant 
grass (Han et al., 2008; Zhou et al., 2009). Quickly establishing more protected areas containing 
S. bungeana grassland is also highly recommended to slow the fast decline in its distribution area. 


5 Conclusions 


The S. bungeana grassland is an invaluable native ecosystem that originally was widely 
distributed on the Loess Plateau. But more than half of its distribution area has been lost over the 
past 50 a due to human activities, rendering it an endangered ecosystem. Quick action is needed 
to restore and conserve this grassland ecosystem through revegetation projects and the 
establishment of protected areas. This modeling study also highlights how an assessment within 
the Red List of Ecosystems conservation framework can be made effectively in a region of China 
having highly fragmented vegetation and a long history of human use. 
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